Impurity Effects on the Nodal Structure of Anisotropic Superconductors 
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For anisotropic superconductors, the gap function belongs to an irreducible representation of the 
point group of crystals, which determines the nodal structure of the gap function. Impurity effects 
on anisotropic superconductors have been treated by the Born or T-matrix approximation by which 
only a partial sum of the pertutbation series is calculated. Instead, we take into acount all terms of 
the perturbation series. As a result, we find that the introduction of nonmagnetic impurities does 
not change the nodal structure in all cases. We also find how to obtain nodes of the gap function 
and show the list for Oh, D4.i1, and D^h- 

PACS numbers: 74.20.Rp, 
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Measurements of direction dependent specific heat and 
direction dependent thermal conductivity has enabled us 
to investigate concrete nodal structures of anisotropic su- 
perconductors. For anisotropic superconductors, the gap 
function is considered to belong to an irreducible rep- 
resentation of the point group of crystals, which deter- 
mines the nodal structure of the gap function-^— . The 
nodal structure of the gap function has strong influence 
on low temperature behaviors of physical quantities, for 
example, specific heat and NMR relaxation rate, which 
arc called the power-laws. Thus, from qualitative point 
of view, it is more important to investigate nodal struc- 
tures than to culculate concrete values. Furthurmorc, 
since what an irreducible representation the gap function 
belongs to has close connection with pairing interaction, 
its identification is crucial for understanding the pairing 
mechanism. 

For conventional s-wave superconductors, low temper- 
ature behaviors are insensitive to small concentrations of 
nonmagnetic impurities (Anderson's Theorem^). How- 
ever, it is not the case with anisotropic superconductors. 
Particularly, for gaps with line nodes, density of states 
(DOS) at the Fermi energy is shifted to nonzero value, 
and therefore low temperature behaviors arc modified by 
arbitrarily small impurity concentrations^. Thus, im- 
purity effects play an important roll to identify the rep- 
resentation of the gap function. 

In this paper, we deal with nonmagnetic impurity ef- 
fects for anisotropic superconductors. We assume that 
the gap function as an order parameter is obtained by the 
generalized Ginzburg-Landau theory (GL). Previously, 
the correction of self-energy part in the Dyson equation 
(Gor'kov equation) due to impurities has been treated by 
the Born or T-matrix approximation. Instead, we take 
into acount all terms of perturbation series of self-energy 
part, not as an approximation. We prove that the intro- 
duction of nonmagnetic impurities does not change the 
representation of the gap function in section IV. We also 
show how to obtain nodes of the gap function and show 
the list for Oh, D§h, and D±h in appendixes. 



We consider a mean field approximated Hamiltonian 
which includes nonmagnetic impurity term, 
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Here, k is a wave number of electrons, s is pseudospin of 
electrons, e(k) is the band energy, V SlS2S3S4 (k, k') is an 
effective electron-electron interaction, V is volume of the 
system, and a ks (a ks ) are annihilation (creation) opera- 
tors. A ss /(fe) is the gap function. The last term is non- 
magnetic impurity term^i By (• ■ ■ ) we mean not only 
expectation value at temperature T, but also disorder- 
averaging. 

Then, equations which temperature Green's functions 
satisfy (Gor'kov equation) are 

G{k,iu> n ) =G (k,iuj n ) 

-G (k, iuj n )±^(k, iu n )G{k, Mn ) (3) 
-G (fc, iu n )£^ (k, iu} n )F{k, -iw n )+, 
F(fc,-w„)t = - G (-fc,iw„)*S (2) (-fc,^ n )*G(fc,iw„) 

-G (-k, iLJ n )*±W{-k, iuj n )*F(k, -uo n )\ 

(4) 

where G(k,i(d n ) and F(k,iuj n ) are normal and anoma- 
lous temperature Green's functions of matrix form re- 
spectively (subscript "0" means that for free fermion sys- 
tem), uj n is the Matsubara frequency, and £W(fc,WJ„) 
and & 2 \k,iu) n ) are the proper self-energy parts. 
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FIG. 1. Gor'kov equation 



comparison with Eq.©, we find that, instead of Aj(fe), 
we can treat ijj(k) or d(k) which satisfy: 



singlet: 



tf(fc) = X>^ (a) (fc), (8) 

i 

4 a) (R-'k) = ^ a) (k)D^(R), (9) 



triplet: 



d(fc)=5>d< a) (fc), 



(10) 



R-'dt\R- 1 k) = Y J df\k)D^(R), (11 



Note that F(k, — iuj n )^ is not merely complex con- 
jugate of F(k,—iLu n ), but hermitian conjugate of 
F(k, —iuj n ) which differs from that in the appendix of 
ref. 8. Throughout this paper, we take an appropriate 
regularization if needed. 



where R is a rotation matrix which does not have the 
parity transformation part of R. Throughout this paper, 
we take the basis functions (k) and d\ a ' (k) to be real 
because accidental degeneracy may not occur. 

Now, the quasiparticle energy depends on A(fe) as 



III. PROPERTIES OF THE GAP FUNCTION 
OF CLEAN SUPERCONDUCTORS 

If temperature T is very close to a transition point 
T c , the gap equation becomes a linear equation. There- 
fore while T is close enough to T c , the gap function is a 
combination of basis functions of an irreducible represen- 
tation r a of the group G (the largest group which does 
not change lattice), such as 



singlet: 



A(fc) = 5>A (Q) (fc)- 



(5) 



A| (k) is i-th basis of irreducible representation T a . 
These coefficients {ci} are determined by GL. Note that 
basis functions are the ones for following transformation 
which forms faithful representation of G: 

U R Lt\R- 1 k)Ul = Y,^\k)Df{R) (6) 

o 

U R = cxpH<7 • n.e/2) (7) 

R is a transformation of G. <r is the Pauli's spin matrix. 
8, n are the rotation angle and the unit vector which rep- 
resent the direction of the rotational part (in the sense of 
right handed screw) of the transformation Rr 1 . dI"\r) 
is an (ij) component of a matrix of irreducible representa- 
tion T Q of G for R. Transformations for pscudospins are 
derived from the fact that due to spin-orbit coupling, in- 
teraction V SlS2SaSi (k, k') is not invariant under rotations 
of wavenumbers unless we rotate pseudospins simultane- 
ously. 

We can write the gap function with ip(k) or d(k) corre- 
sponding to singlet (even parity) or triplet (odd parity), 
such as A(fc) = i(Tytp(k), A(k) = i(d(k) ■ cr)a- y , respec- 
tively, where <x and a y are the Pauli's spin matrices. By 



E k = ;/e(fc) 2 + -trA(fc)A(fc)t = ^e{kf + \4>{k)\\ 

(12) 



triplet: 



E k± = y/e(k)* + \d(k)\*±\q(k)\, 



(13) 



where q(k) = id(k) x d(k)* . A(fe) has zero points de- 
rived from its belonging irreducible representation (See 
Appendix B). The zero points make nodes of the excita- 
tion energy on the Fermi surface. 



IV. SUPERCONDUCTORS IN THE PRESENCE 
OF NONMAGNETIC IMPURITIES 



Since we can find that the quantity in the presence of 
nonmagnetic impurities corresponding to the gap func- 
tion of clean superconductors is the proper self-energy 
part of anomalous type S' 2 )(fc, iw n ), let us consider 
£( 2 )(fc, iuj n ). For convenience, we first consider a term 
GxFiG* of perturbation series of S( 2 )(fc,iw„) (SeeFigHJ. 
Note that the subscript "1" means the exact Green's 
functions in the absence of impurities. 
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FIG. 2. 
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The transformation of T,^ 2 \k,ico n ) such as 
£W(k,u) -> UatMtRr^iUrjUZ leads to corre- 
sponding transformation of G\F\G\ term included in 
& 2 '(k,iw n ) such as 

n 2 J2 V(fc-fei) 2 V(*i -k 2 ) 2 

fel ,fc 2 

x Gi{k -fei + k2,iuj n )Fi(k 2 ,iuJn)G* 1 (-ki 1 iuj n ) 

(14) 

Yl V(R' 1 k-ki) 2 V(k 1 -k 2 ) 2 
fci ,fc 2 

x U R G 1 (R~ 1 k - fci + k 2 ,iuj n ) 
x A(fc,to n )GI(-fci,iw n )E)£ (15) 
=n? ^ F(ir x fc - ie _1 Jfei) a V'(JJ- 1 fci - iT 1 *^) 2 

x U R G 1 (R- 1 k - i? _1 fei + ir 1 fc 2 ,iw„) 

x A(^" 1 fe2,iw„)GJ(- J R- 1 fe 1 ,^„)^ (16) 

=n 2 ^ ^(iJ-^fe-fei)) 2 ^- 1 ^! -fc 2 )) 2 
fci ,fc 2 

x ^G 1 (i?- 1 (fc-fc 1 +fc 2 ),^„)[/^ 
x U R F 1 (R- 1 k 2 ,icJ n )UR 

xUZG* 1 (-ir 1 ki,iu n )UZ, (17) 
now, if V(fc) satisfies V(iJ- 1 fc) = V(fc), (18) 

=n 2 ^ F(fc-fc 1 ) 2 T/(fe 1 -fc 2 ) 2 

x ^^(^(fe-fei + fea),^)^ 
x U R F 1 (R- 1 k 2 ,icJ n )UR 

xU^GK-R^k^ioj^U^, (19) 

where n., is the impurity concentration. This 
implies that transformation such as S*- 2 -* (fe, iui n ) — > 
UrYj^ (R~ 1 k, iui n )U R is equivalent to the transforma- 
tion of Green's functions in E^ 2 '(fc, iu> n ) such as 

Gi(fc,iw„) -> UrG^R- 1 ^^)^, (20) 

F 1 (k,ioj n ) -> U R F 1 (R- 1 k,ioj n )U^, (21) 

(2?(-fc,iw n ) -> U R Gl{-R- l k,iu n )Ul, (22) 

Ft(fc, _ia, n ) U* R Fl{R- l k, -ioj n )U R . (23) 

Although the current example does not have 
Pl(k,—iuj n ), we can always treat general terms 
like this. The reason for this is as follows. 

Combinations of adjoining two Green's functions ap- 
pearing in S( 2 )(fc,icj„) are restricted as left raw below 
because of particle number conservation at each point. 
And we can insert spin rotation matrices into two Green's 
functions as right raw below. Moreover, it is always 



possible to transform internal wave numbers as well as 
G\F\G\ because internal wave numbers are merely inte- 
gral variables . 

FxG* -> F 1 U T U*G* V 
FiP} -> F 1 U T U*Pl 
->■ GIU T U*G{, 
G{Pl -> &iu t u*f}, 

F+Gi -> P^UGx, 

PIP 1 -> PIU^UFx. 

Therefore examination on what irreducible repre- 
sentations these Green's functions have for above 
transformations (j2T))) - (|2"3"|) will make it clear that what 
irreducible representations E' 2 ' (fe, io; n ) has for above 
transformations for & 2 ^ (k,ioj n ). 

Here, we define G', T'g, and D^. By G' we mean the 
maximal subgroup of G where the quasiparticle energy 
of clean superconductor is invariant. We can find that G' 
is not only maximal but also the largest. Because if G" 
where the excitation energy is invariant is not contained 
by G', the quasiparticle energy is also invariant under 
the subgroup which is generated from G' U G"(D G'), 
this contradicts the fact that G' is maximal. L^ is an 
irreducible representation subspace of G'. 

D'p is a subspace of T'p generated from elements which 

have a form of ^\ c^tp'^ = Yli^^l ■ The coefficients 
{ci} are determined by GL. {c^} are determined uniquely 
from {ci}. 

Let us examine what irreducible representations 
Green's functions have under above assumption. In what 
follows, A(fe) is called unitary if the product A(fe)A(fe)' 
is proportional to the unit matrix, otherwise it is called 
nonunitary. Note that only triplet case can be nonuni- 
tary. 

It can be proved that regardless of the absence or the 
presence of impurities, A(fc) € D'^ C T'a until level cross- 
ing (including by another bifurcation from the solution) 
occurs for any GL solution^. In what follows, we shall 
consider two cases (i) and (ii) separately. 

(i) In case order parameter coefficients {ci} obtained as 
a solution of GL are all real (complex number multiplied 
overall is neglected), only unitary case applies to this 
case. 

Exact Green's functions of clean superconductors are 

Gi(k,iw n ) = 9 a , (24) 

< + E k 

fi l (k,i Un ) = --£&L. (25) 

Because the excitation energy of clean superconductors 
and the band energy satisfies = e(fe') = e(fe) 



4 



under transformations of G respectively, Green's func- 
tions are transformed by G' as 



U R F 1 {R- 1 k,iw n )Ul 



G 1 (k,iuj n ), 
U R K{R- l k)tJl 



(iR - lk) -URV U R 



Hence, 

Gi(fe,™ n )er'+ 

6 D'p. 

Because {ci} are all real, 

[Dp <8> D^Jxs = 0, 



(26) 

(27) 

(28) 
(29) 

(30) 



where [• • • ]as means antisymmetric tensor product. Fur- 
thermore, since we can find that the condition which the 
excitation energy for unitary case is invariant under G' 
is equivalent to 



Therefore, 

D'o(g 



[D'^D'^scT't 



D' f ; = [D' fj ® D>f] s + [D'p 

cr' 1 + + o = r'+. 



Dp] AS 



(31) 



(32) 



Now, the representation of each term of Y,^ 2 \k, iuj n ) is 




Dp c D'p. 



(33) 



Here, we neglect G\{k,iw n ) and G\(k,iw n ) since they 
belong to r^ + . 

And the facts that F(k,iw n ) and F'(k,— iuj n ) appear 
in turn, and the number of F(k,ioj n ) is always greater 
than F^k, —iuj n ) by one are used. 

Because all terms appearing in E( 2 )(fc,iw„) belong to 
the same D'p as the original A(fe), if equation (|32|) holds, 
E*- 2 -* (k, iu> n ) belongs to the same D'p as the original A(fe) 
under G' . It should be noted that in the present case, 
in the table in Sigrist-Ueda£, T'p are all ID irreducible 
representation, hence Dp = T'p. 

(ii) Whereas, in the case of complex coefficients {cj}, 
equation (|3"2"j) does not hold. Because there is no such ID 
representation of G, it is enough for us to prove the case 
of 2D and 3D. 

For G' , any gap function which has complex coefficients 
is 2D irreducible representation^. Now, [D'p® D'p*]as is a 



certain ID representation other than because of its 
dimension and it can be proved that it is not T'^~ . Hence, 
we find 

[D'p®Dp*]AsCT' 1D) (34) 

where T' 1D is defined as a certain ID representation other 
than r'^ which has [Dp ® Dp] as as a subset. 

It is worth mentioning that in this case, nonunitary 
cases arise, and it is easy to verify that 

£i(fc,iwn) G r'+ + [D'p ® D'p*] AS , (35) 
Fi(k, iu> n ) e D'p, (36) 
in nonunitary cases. 

In any case of (ii), i.e., unitary or nonunitary, general 
term of E^ 2 ^(fe, iu> n ) has the following representation: 



D'p + T' 1D ® D'p 



(37) 
(38) 



We can find that in all cases, this is Dp, and the solution 
of the gap equation is not affected by the introduction 
of nonmagnetic impurities^. 



After all, in any case, we can find that £( 2 )(fe, iui n ) G 



D' R 



V. SUMMARY AND DISCUSSION 

We proved that the introduction of nonmagnetic impu- 
rities does not change the nodal structure of anisotropic 
superconductors. We also show why zeros of S*- 2 ^(fc,0) 
are related to the nodal structure in Appendix A, how to 
obtain nodal structures in Appendix B and show them in 
the case of O^, D 6 h, and D±h in the table II, III, and IV 
of Appendix C respectively. Calculations in Appendix C 
and D were done by using MATHEMATICA. 

Although theoretical culculations for anisotropic s- 
wave gaps show that the gap anisotropy is smeared out by 
introducing nonmagnetic impurities^i, and anisotropic 
s-wave superconductivity is therefore anticipated to ex- 
plain experimental evidences of YW^^G^—, it is no 
wonder that the nodes disappear because the nodes are 
not the ones as an irreducible representation. 

We should pay attention to that even if nodes of 
gaps disappear due to impurities, we do not say any- 
thing about the amplitudes of gaps. Hence, the fact ob- 
tained in this paper does not contradict the fact that the 
axial-like gap (gap with point node(s)) and the polar- 
like gap (gap with line node(s)) are suppressed partially 
by introducing nonmagnetic impurities in the case of the 
Born approximation^, and in the case of the T-matrix 



approximatioi 



.15-17 



. Especially, a theoretical culculation 



by Haas et alr^ and an experiment by Hashimoto et alJ£ 
agrees with our theory, although these are too simple to 
be evidences. 
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Appendix A: The Existence of the nodes 

From Gor'kov equation, we can derive that 

((*w„ - e(k))& - ± (1) (k,iLu n ) + ^(k^lLJn) 

X ((-iw„ - c(fe))<7 + liOn)*)- 1 

x £ (2) (-fc,iu; n )*)G(fc, 

= -to (Al) 

Putting ui to 0, we can obtain zero energy excitation. 
Consider a point where £( 2 )(fc,0) = 0. At that point, if 
((— e(fc))cf + X^ 1 ^— fc, 0)*) -1 is divergent more strongly 
than two £( 2 )(fc,0)s on both side, the Green's function 
has no pole there. In this case, there is no node. 

If it is divergent less strongly than or equal to two 
£( 2 )(fc, 0)s on both side, or not at all divergent, the equa- 
tion becomes 



F 1 (k 2 ,0)G 1 {-k 1 ,0) T 
= F 1 (k 2 ,0)G 1 (-k 1 ,0)*, 

A(fei,0)t(-G 1 (fc 2 ,0))^ 
(-Gi(fe2,0))Fi(fci,0) = 
(-Gi(fc2,0))Fi(fci,0), 

A (fe 1; 0)tFi(fc 2 ,0)-> 
F 1 (k 2 ,0yP 1 (k 1 ,0). 

Note that in E*- 1 -* (fe, 0), the number of F\ and F{ is 
cquall. Similar to the argument of E^ 2 ^(fc,0), the trans- 
formation causes a change with dually corresponding 
term with no sign change. 

Thus, EW(fe,0)t = SW(fe,0), and (e(fe)a + 
£W(fe, 0)) can be diagonalized by a unitaty matrix. 

Let us write the eigenvalues as e(fe) + E^(fe,0), (i = 
1,2). 

e(fe) changes its sign on the Fermi surface, while 

£^ (fe, 0) should be slowly varying in the vicinity of 
the Fermi surface if impurity concentration is sufficiently 
small although this needs to be proved. Therefore, 
there should be a point near the Fermi surface where 
e(fc) + £-^(fc, 0) = 0. At that point, zero energy excita- 
tion exists. 



((e(fe)<7 +XI (1) (fe,0))G(fe,0) 
= -<5o, 



(A2) 



where we renormalize the residue into EW(fe,0) if it 
exists. 

Moreover, EW(fe,0)t = EW(fe,0) because the 
transformation £W(fc,0) -> £( 1 >(fc,0) t causes transfor- 
mations on internal Green's functions such as, 

(-Gi(fei, 0))(-Gi(fe 2 ,0))-». 
(-G 1 (fe 2 ,0))t(-G 1 (fe 1 ,0))t 
= (-Gi(fc 2 ,0))(-Gi(fci,0)), 

(-<5i(fci,o))A(fc2,o)-> 

A(fc 2 ,0)t(-G 1 (fc 1 ,0))t 
= F 1 (fe 2 ,0)(-G 1 (-fei,0)), 

F 1 (k 1 ,0)G 1 (-k 2 ,0)* ->■ 
G 1 (-fe 2 ,0)*F 1 (fei,0)t 
= G 1 (-fe 2 ,0)*Fi(fei,0)t, 

A(fci,o)A(fc 2 ,o)t^ 
A(*2,o)A(fci,o)t, 

G 1 (-fc 1 ,0)*G 1 (-fc 2 ,0)* -> 
G 1 (-fe 2 ,0) T Gi(-fei,0) T 
= G 1 (-fe 2 ,0))*G 1 (fe 2 ,0)*, 

G 1 (-fe 1 ,0)*F 1 (fe 2 ,0)t^. 



Appendix B: How to obtain nodes of the gap 
functions 

The fact that £( 2 )(fc,0) has the same parity as the 
gap function is almost trivial. This condition and the 
condition £( 2 )(-fc,0) T = -£( 2 )(fc,0) leads to the fact 
that E' 2 '(fe,0) can be written in the same form as the 
gap function, i.e., E( 2 )(fc,0) = ia y %j){k) or £( 2 )(fc,0) = 
i(d(k) ■ <r)(Jy with certain functions ip(k) or d(k), corre- 
sponding to the parity, respectively. 

Searching for zeros of £( 2 )(fc,0) is reduced to 
searching for those of ip(k) or d(k). The fact 
£( 2 )(-fc,0) T = -£( 2 >(fc,0) is proved in the latter 
half of this appendix. 

Although there are many sets of basis functions of 
irreducible representations, it may occur that any basis 
functions become zero at certain points. Here, we 
consider such zero points. 

singlet case: 

We write basis functions {Vv }i=i,2,— ,d a of the irre- 
ducible representation T a (d a dimension) of the point 
group G in the form: 



^(fe) 



4 a) 



(Bl) 
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And ip(k) satisfies 

ipiR^k) = D{R) T ip(k), 



(B2) 



where (D(R))ij = D^\r) is a representation matrix of 
the group element R of the irreducible representation T a . 

The equation (|B2|) gives relationships between values 
of ip(k) at points which is moved by elements of G. 

Now, let us consider a point k which is not invariant 
under any element R of G. Then, {i? _1 fc}^ 6 G consists 
of points of the same number g as the order of G, and 
values of ip{k) at {i? _1 fe}ij e G have relationships. Since 
there are at most d a (g — 1) independent linear equations 
for d a g variables {t/>(.R _1 fc)}ii e Gas 



n 



i 



-D(R2) T \ 
-D(R 3 ) T 



i -D(R g r; 

the solutions are 



f^(R^k)\ 
ipiR^k 

^(Rg'k 

v m j 



= 0, (B3) 



^(R^k) 



v m J 



f(D{R 2 ) T \ \ 

?, 



D(R^ T 



D(R 9 ) T 

V i J 



a 



(B4) 
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where a is an arbitrary d a dimensional vector, and 
D(R\) — 1 (d a dimensional unit matrix). Therefore no 
restriction that forces {■0(i? _1 fc)}ij 6 G to be zeros ap- 
pears. 

Next we consider a point k which is invariant under 
some elements of G. These points can exist only when 
they are on rotational axes or mirror reflection planes. In 
contrast to previous case, since ip(k) with the condition 
R~ x k = k satisfies 



ip{k) = D(R) T i>(k) 
(D(Rf - i)V(fc) = 0, 



(B5) 
(B6) 



the value of ip(k) is restricted. 

(i) In the case that there exist R such that the rank of 
the matrix (D(R) T — 1) equals d a , that is, determinant 
of the matrix is not equal to zero, then tp(k) = 0. In 
this case, for any set of coefficients {cf}, ^ZjCfV'f (k) — 
c-V(fe) = 0. For example, (0,0,1) (line) and [0,0,1] 
(point) of r 5+ of D^h applies to this case. 

(ii) In the case that there is no R such that the 
rank of the matrix (D(R) T — 1) equals d a . In this 
case, if the rank of the matrix for an R equals n < d a , 
the solutions of the equation (|336|) must lie in d a — n 
dimensional subspace. If k has several R which satisies 
R~ 1 k = k, the subspace of solutions of the equation is 
the common subspace of subspaces for each R. i/)(k) 
is not always in this case. However, if coefficients 



c = (c"c2 ■ • ■ c°[ ) T belongs to the orthocomplcment of 
the common subspace, then c • ip(k) = J2i c ? V'f (^) = 0- 
This means that the gap function is zero at the point k. 

triplet case: 

The same manner as singlet applies to triplet case ex- 
cept for little complexity. Basis functions {df }i=i,2,--- ,d a 
of the irreducible representation T a satisfies 



R-^t\R- 1 k) = Y J df\k)D^{R) 



(B7) 



where R is that which parity transformation is reduced 
from R. We write 



( d i ] \ 



v 4:V 



(B8) 



In this notation, at the point k with the condition 
R~ 1 k = k, we have 



(R- 



R- 



'Du(R) D 2 i{R) ■ 
d(k) = I D 12 (R) D 22 (R) I d(fe ) 




D^R) D 21 {R) 
D 12 (R) D 22 (R) 



(B9) 

> d(k) = 0, 
(BIO) 



I 0\ 

where D l} (R) = D\f(R) 1 and I is a 34 di- 

V° ° V 

mensional unit marix. In the case corresponding to (ii) 
of singlet, we must note that 



/ d^\ 














0- ■ 


r a 







4 Q) 


r? 
' 1 








' 2 


0- • 





r° 
c d a 







r" 
1 1 








c 2 


• 





c dJ 


1 4:v 



= ^2c?d^(k). 

(Bll) 

Since discussion on triplet case will be the same as 
singlet case except for this, we abbreviate furthur detail. 

The proof for £ (2 >(-fc,0) T = -£ (2) (fc, 0): 

A transformation £ (2) (fc,0) -> S (2) (-fc,0) T induces 
transformations on internal Green's functions such as: 
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(-G 1 (fc 1 ,0))(-G 1 (fc 2 ,0))^ 
(-G 1 (-fe 2 ,0)) T (-G 1 (-fc 1 ,0)) T 
= G 1 (-fe 2 ,0)*G 1 (-fe 1 ,0)*, 

(-Gi(fci,0))Fi(fc2,0)-> 
F 1 (-k 2 ,0) T (-G 1 (-k 1 ,0)) T 
= -F 1 (k 2> Q)G 1 (-k 1 ,0)*, 

A(fci,0)G 1 (-fc 2 ,0)* ->■ 

Gi(fe 2) 0)tF 1 (-fei,0) T 

= G 1 (-fe 2) 0)*(-l)A(fei,0), 

A(fei,o)A(fca,o)t -> 

Fi(-fc 2) 0)*Fi(-fci,0) T 
= F 1 (fe 2) 0)tF 1 (fei,0), 

G 1 (-fe 1 ,0)*G 1 (-fc 2 ,0)* -> 

Gi(fc 2 ,0)tGi(fc 1 ,0)t 

= (-Gi(fc 2 ,0))(-Gi(fc 2 ,0)), 

G 1 (-fc 1 ,0)*A(fc 2 ,0)t^ 

A(-fe 2 ,0)*Gi(fe 1 ,0)t 

= (-l)A(-fe 2! 0)tGi(fei,0), 

A(fci,0)t(-G 1 (fc 2 ,0))^ 
(-G 1 (-fe 2 ,0)) T F 1 (-fci,0)* = 
(-Gi(fea,0))(-l)Fi(fci,0)t, 

Fi(fei,0)tFi(fe2,0) -)■ 

A(-fc 2 ,o) T A(-fci,o)* = 

A(fc 2 ,0)Fi(fe 1 ,0)t. 

Here, the transformation of external wave number k — > 
—A; is equivalent to the transformation of internal wave 
number ki — > — fej. 

Thus, by the transformation S^(fe, 0) — )• 
S (2) (-fc,0) T , each term in £( 2 )(fc,0) is multiplied 
by a factor —1, because —1 is multiplied "the summa- 
tion of the number of F\ and P\ " times, and they are 
always odd. 

In the case of a diagrammatically symmetric term, the 
transformation causes only the change of internal wave 
number indices, and thus —1 is multiplied. 

In the case of a diagrammatically asymmetric term, the 
transformation causes a change with dually correspond- 
ing term, and thus —1 is multiplied. 

Thus, £( 2 )(-fc,0) T = -£ (2) (M). 



Appendix C: Table of Nodal Structures of Oh, Dqh, 
and Anisotropic Superconductors 

We choose basis functions of irreducible representa- 
tions for Oh,Deh,aiid D^h as follows. 



1 I: basis functions of irreducible representations for Oh, D^h, and D^h 

r V'(fc)/d(fc) 



o h 

parity even (singlet) 







(fe) 


= 1 


1 2 




(fe) 


= (fe 2 -fc 2 )(fe 2 -fc 2 )(fe 2 - fe 2 ) 
V x y I \ y z I V z x > 


It 




(fe) 


= 2fc " - k l - k l 






(fe) 


= WZ-k'v) 




< +) 


(fe) 


— k y k z (k y — k z ) 






(fe) 


— k z k x (k z — k x ) 




<4 4+) 


(fe) 


— kxk y (k x — k y ) 


rt 




(fe) 






^ 2 5+) 


(fe) 


= fclfel 






(fe) 


= k x ky 

parity odd (triplet) 
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(fe) 


— xk x -\- yky -\- zk z 




d (2 -' 


(fe) 


= xfc*(fc 2 - k 2 y )+yk y + (fe 2 - fe 2 ) +zfc,(fc 2 - fe 2 ) 


r~ 

3 


d< 3 -» 


(fe) 


— 2zk - xk - yk 




d< 3 -' 


(fe) 


— \/~3(xk x — yky'j 


1 i 




(fe) 


— nk — v k 






(fe) 


— zk x — xk ~ 






(fe) 


— Ojt\,y yf^x 


1 5 




(fe) 


— yk z ~\- zky 






(fe) 








(fe) 


— xk,, -\- lifc™ 

tXjn.y 7 y^X 








D 6h 








parity even (singlet) 






(fe) 


= 1 


r+ 


^(2 + ) 


(fe) 


= fexfe„(fe 2 - 3fe 2 )(fe 2 - 3fc 2 ) 




^(3 + ) 


(fe) 


— kzk x (k x 3fcy) 


r 4 + 




(fe) 


— k z ky{ky 


r 5 


i>f + ) 


(fe) 


_ k ^ kr 




4 5+) 


(fe) 


= kykl 


r+ 




(fe) 








(fe) 






d (1 -' 




parity odd (triplet) 


r _ 
1 1 


(fe) 


— t 1 fe* —I— itk 1 ^ k* 


r _ 


d< 2 -' 


(fe) 


— A* y — y A.' 


A 3 


d< 3 -' 


(fe) 


= zfc (fe 2 — 3fc 2 ) fe (fe 2 — fe 2 )x — 2fc fc -u 


rr 
i 




(fe) 


— zk y (ky - 3k^), k z (k y — k^)y — 2k x k v x 


1 5 


j(B-) 
a i 


' fe' ) 


= xk z ,zk x 




j(B-) 

a 2 


' fe) 


— yk z , jsfcy 


r~ 


a. 


(fe) 


— T* fe 1 — i i k' 

— tL, " J X, y "- y 






(fe) 


— xky -f- yk x 
















parity even (singlet) 


r+ 
1 1 


,/,(!+) 
■r 


(fe) 


— i 


r+ 


^(2 + ) 


(fe) 


= k-r-ki.lk 2 - k 2 ) 

n, x n, y\ n 'x yJ 


r+ 


v,(3 + ) 


(fe) 


— k 2 — k 2 

^x ^y 


rt 

1 


^(4 + ) 


(fe) 








(fe) 


= k x k z 






(fe) 










parity odd (triplet) 




(fe) 


— xk x + yky , zk z 




d< 2 -> 


(fe) 






d' 3 -' 


(fe) 




rr 


d< 4 -' 


(fe) 


— xky -\- yk x 




df"> 


(fe) 


— xk z , zk x 




df > 


(fe) 


— yk z , zky 



Order parameters for Oh, Dqh,, and D^h obtained from 
GL and the change of their nodal structure following 
the introduction of nonmagnetic impurities are given 
in the following tables. For example, (l,w,u; 2 ) is a 
combination of expansion coefficients for basis functions 
defined here, where u> = e' 27r / 3 . (a) means numbers of 
order parameters being solutions of GL simultaneously, 
(b) means the largest subgroup G' where the excitation 
energy is invariant (Ty ). (c) means zeros in fc-space 
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immediately under the transition point T c in the absence 
of impurities. [k x , k y , k z ] in the table means a straight 
line generated by a vector (k Xl k y ,k z ). (k x , k y ,k z ) in the 
table means a plane which is perpendicular to a vector 
(k x ,k v ,k z ) and goes through the origin of fc-space. (d) 
means zeros below T c in the absence and the presence of 



impurities. Nodes are generated where zeros and Fermi 
surface intersect, "none" means nonexistence of node, 
and left arrow means that the nodes on the left-hand side 
does not disappear. In the right edge, which representa- 
tion appears is written in the case of complex coefficients. 



TABLE II: O h 



(o h )r 


</i(fe)/d(fe) 


(a) 


(b) 


(<0 


(d) 


rt 


1 


1 


o h (r+) 


none 


none 


it 


1 


1 


o h (r+) 


(1,±1,0),(0,1,±1), (±1,0,1) 




rt« 


(0,1) 


3 


p,(0,0, 1) \ 

-^4/1 V 1 3 J 


(1,±1,0) 




r+(2) 


(1,0) 


3 


n (0,0.1) ,„+., 


[1,1,1], [-1,1,1], [1,-1,1], [1,1,-1] 


none 


r+( 3 ) 


CM) 


2 


o h (r+) 


[1,1,1], [-1,1,1], [1,-1,1], [1,1,-1] 


<~ 


rj(i) 


(l,w,w 2 ) 


8 


D< 1 d 1 ' 1 '(r+) 


[1,1,1] 




rt(2) 


(1,1,1) 


4 




(1, -1,0), (0, 1, -1),(-1,0, 1) 
[1,-1,0], [0,1,-1], [-1,0,1] 


<- 
*~ 


rj(3) 


(1,0,0) 


3 


T-)(1,0,0) /p + \ 
^4h V 1 2 J 


(0,1,±1),(0, 0,1), (0,1,0) 


*~ 


It (4) 


(l,i,0) 


6 


p. (0,0, 1) /p+ * 
^4h V 1 5 J 


(0,0, 1),[0,0, 1] 


*~ 


it« 


(l,w,w 2 ) 


8 


< 1 ' 1) (r+) 


[1,1,1] 

[1,0,0], [0,1,0], [0,0,1] 


<- 


it (2) 


(1,1,1) 


4 


< ia, (r-) 


[1,0,0], [0,1,0], [0,0,1] 


none 


it (3) 


(1,0,0) 


3 


p. (0,0, 1) \ 

'-'ih 4 ; 


(0,1,0), (0,0,1) 




it (4) 


(l,i,0) 


6 


n (0,0.1) ,p+x 
^4h V 1 5 J 


(0,0, 1),[0,0, 1] 


<- 


17 


1 


1 


Oh (IT) 


none 


none 




1 


1 


Oh (17) 


[1,0,0], [0,1,0], [0,0,1] 




r 3 "(i) 


(0,1) 


3 


DS 0,1) (r 3 -) 


[0,0,1] 




^ (2) 


(1,0) 


3 


\ L 1 I 


none 


none 


r 3 (3) 


CM) 


2 


Oh(r 3 -) 


none 


none 


17 CD 


(l,^,^ 2 ) 


8 


< 1 ' 1, (r 3 ) 


none 


none 


17(2) 


(i,i,i) 


4 


^ 1 /' 1) (r 2 -) 


none 


none 


rr(3) 


(1,0,0) 


3 


^4h V 1 2 1 


none 


none 


ri(4) 


(0,-t,l) 


6 


D a.o.o) (r - ) 


none 


none 


r 5 "(i) 


{l,u> 2 ,u>) 


8 


^ 1 /' 1, (r 3 -) 


none 


none 


r 5 "(2) 


(1,1,1) 


4 




none 


none 


r 5 " (3) 


(1,0,0) 


3 


D (i,o,o), r -, 

^4h V 1 4 -> 


[1,0,0] 




I7W 


(0,i,l) 


6 


^4h \ L 5 J 


none 


none 



TABLE III: D, 



(D 6h )r 


i>(k)/d(k) 


(a) 


(b) 


(c) 


(d) 


it 


i 


1 


D 6h (T+) 


none 


none 
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TABLE III: D 6h 



(D 6h )V 


V-(fe)/d(fc) 


(a) 


(b) 


(c) 


(d) 


r+ 

1 2 


1 


1 


-D6/i(r^) 


(1, 0, 0),(0, 1, 0),(1. ±\/3, 0),(±\/3, 1, 0) 




r i 


1 


1 


-^6h (^3") 


(1,0,0), (0,0,1), (1,±V3,0) 






1 


1 


Deh(Pt) 


(O.l.OJXO.O.lJ^iv^, 1,0) 




r 5 + (i) 


(1,0) 


3 


^2fe(r+) 


(1,0,0), (0,0,1) 




r 5 + (2) 


(0,1) 


3 


^2h(r+) 


(0,1,0), (0,0,1) 




r+(3) 


CM) 


2 


•Deh(r^) 


(0,0,1), [0,0,1] 




r tW 


(1,0) 


3 


£>2/.(r+) 


[0,0, 1] 


none 


r?(2) 


(0,1) 


3 


-D2/.(r+) 


(1,0, 0), (0,1,0) 




rJ-(3) 


CM) 


2 


^6/i (^6 ) 


[0,0,1] 
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1 


1 


D 6h (r-) 


none 




r~ 


1 


1 


Den (1^2 ) 


none 






1 


1 


^6h (To ) 


[0,0,1] 




r I 


1 


1 


D 6h (r7) 


[0,0,1] 




rj(i) 


(1,0) 


3 




none 


none 


r 5 "(2) 


(0.1) 


3 


A^OT) 


none 


none 


Ts (3) 


CM) 


2 


A5h(r 5 ~) 


none 


none 


r 6 -(i) 


(1,0) 


3 


^2h(rr) 


[0,0,1] 


none 


r 6 "(2) 


(0,1) 


3 


^2h(r-) 


[0,0,1] 


none 


Te (3) 


CM) 


2 


-D6/i(r 6 ) 


[0,0,1] 





TABLE IV: D ih 



(D ih )T 


V(fe)Mfe) 


(a) 


(b) 


(c) 


(d) 


rt 


1 


1 


£> 4fe (r+) 


none 


none 


r 2 + 


1 


1 




(1,0,0),(0,1,0),(1,±1.0) 


•J— 


?t 


1 


1 


o 4 h(r+) 


(1,±1,0) 




rt 


1 


1 


£Mr+) 


(1,0,0), (0,1,0) 




r 5 + (i) 




2 


£>4h(r+) 


[0,0,1], (0,0,1) 




r 5 + (2) 


(1,1) 


2 


n (diagonal) , p +v 
W 2k ^ 3 J 


(0,0,1), (1,1,0) 




r 5 + 0) 


(1,0) 


2 




(1,0,0), (0,0,1) 




rr 


1 


1 




none 


none 


r~ 


1 


1 




none 


none 


r 3 " 


1 


1 




[0,0,1] 




ri 


1 


1 




[0,0,1] 




i7(i) 


CM) 


2 




none 


none 


r 5 "(2) 


(1,1) 


2 


D (diagonal) g,-) 


none 


none 


r 5 "(3) 


(1,0) 


2 


C2h(r 3 -) 


none 


none 



Here, (0,0,1) of D^' u ' 1J (r 3+ ) means as follows. O h 
has three kinds of subgroup D^h which have different four 
times rotation axis each other, for example, D^f' 1 ^ (T^) 
is a D^h which has (0, 0, 1) as four times rotation axis. 
Although such subgroups are different for every order 
parameter, we write only subgroups for representative 



order parameters. 

The largest subgroup where the quasiparticle energy of 
a given order parameter is invariant can be obtained as 
follows^ .. First, examine all maximal subgroups whether 
the excitation energy is invariant or not under the 
subgroups. If there exits, that is the largest subgroup. 
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If not, examine all second maximal subgroups, and so 
on. The maximal subgroup obtained like this is also the 
largest subgroup as mentioned in section IV. 



Appendix D: Expansion of Order Parameters of 
Complex Coefficients with respect to G' where the 
Quasiparticle Energy of Clean Superconductor is 
Invariant 

We examined {c' t } of c^f ] = £ 4 c^ (q) . 
means that it belongs to ith basis of /3th irreducible rep- 
resentation of G' where the quasiparticle energy of clean 
superconductor is invariant. These results can be under- 
stood from the GL theory. 



O h 



(2D)cven 

T+ (3) (subgroup G' = O h ) 



r^"(l) (subgroup G' = D^d which has (1,1,1) of Oh as 
a three times rotation axis) 

kyk z — (— Ojk z k x ~\~ k x ky 



[2k x ky k x k z kyk z ) 



- i 3 J 1 3 (v / 3(fc a; - k y )k z ) 

=v ,f+) +# ;(3+) 



(D5) 



(4) (subgroup G' = D ih which has (0,0,1) of O h as 
a four times rotation axis) 



kyk z + ik z k 



(D6) 



2 1.2 



k x + uiky + u> k 



-~ 2 (2kl 



k 2 y ) + kV3(kl-k 2 y )) 



(2D)odd 

^(3) (subgroup G' = O h ) 

xk x + Luyk y + uj 2 zk z = ip'-^ 3 ' + iip'^ 3 ^ 



(Dl) 



(D2) 



(3D)odd 

r^" (1) (subgroup G' = D 3d which has (1,1,1) of Oh as 
a three times rotation axis) 

x(u}k y — u} 2 k z ) + y(k z — uik x ) + z(oj 2 k x — k y ) 
3 + iV3 



4 

-1 + iVS 



(x(-k z ) + y(-k z ) + z(k x + k y ))+ 
(x(-2k v - k z ) + y(2k x + k z ) + z(k x - k y )) 



(D7) 



(3D) even 

r^l) (subgroup G' = Z?3d which has (1,1,1) of Oh as 
a three times rotation axis) 

k y k z {ky k z ) + u)k z k x {k z k x ) + to k x k y (k x fc, y ) 



3-iV3 



(kx ~\~ ky)k z (k x k x ky + k k z )-\- 



(k^i'Zky + k z ) — kyk z (k 2 — k 2 ) — k x (2ky + k z )) 



= ^(3+) +# f+) 



(D3) 



(4) (subgroup G' = D ih which has (0,0,1) of O h as 
a four times rotation axis) 



kyk z (k 2 — k 2 ) + ik z k x (k 2 — k 2 ) 



=v4 (5+) 



(D4) 



F 4 (4)(subgroup G" = D 4 h which has (1,0,0) of Oh as 
a four times rotation axis) 



x(k y + ik z ) + y{~k x ) + z(-ik x ) 



(D8) 



T 5 (1) (subgroup G' = Dsd which has (1,1,1) of Oh as 
a three times rotation axis) 

x(ujk y + Lo 2 k z ) + y(k z + u>k x ) + z(ui 2 k x + k y ) 

3 \x(k z - 2k y ) + y(k z - 2k x ) + z(k x + k y ))+ 



4 

V3 + i 



(V3(xk z + y{-k z ) + z(k x - k y ))) 



(D9) 
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T 5 (4) (subgroup G' = D ih which has (0,0,1) of Oh as 
four times rotation axis) 



x(k y + ik z ) + yk x + z(ik x ) 



(D10) 



IV (3) (subgroup G' = D 6h ) 

(x + iy)(k x + iky) = ' + iip, 



(D14) 



D 



(2D)even 

T+ (3) (subgroup G' = D 6h ) 



k z {k x + ik y )=i)f +) +#| 5+) 



(Dll) 



D 



4/i 



(2D)even 

T+(l) (subgroup G' = D 4h ) 



k z {k x + iky) = ip\ +iip. 



( 5 +) _l_ ,-„/,( 5 +) 



(D15) 



T+ (3) (subgroup G' = D 6h ) 



(k x +ik y f = < +) +z^ 6+) 



(D12) (2D )odd 

T 5 -(l) (subgroup G' = D 4h ) 



(2D)odd 

T- (3) (subgroup G' = D 6h ) 



xk z + iyk z = V4 5 } + i^2~ ] (D16) 



xk z + iyk z = tp[ 5 •* + iipf ^ 



(D13) 
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